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Abstract. We introduce a continuous damage fiber bundle model and compare its behavior with that of 
dry fiber bundles. Several interesting constitutive behaviors, such as plasticity, are found in this model 
depending on the value of the damage parameter and on the form of the disorder distribution. We compare 
the constitutive behavior of global load transfer models, obtained analytically, with local load transfer 
models numerical simulations. The evolution of the damage is studied analyzing the cluster statistics for 
dry and continous damage fiber bundles. Finally, it is shown that quenched random thresholds enhance 
damage localization. 

PACS. 4 6.50.+a, 62.20.Fe, 62.20.Mk 
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1 introduction 

The rupture of disordered media has recently attracted 
much technological and industrial interest and has been 
widely studied in statistical physics. It has been suggested 
by several authors that the failure of a disordered material 
subjected to an increasing external load shares many fea- 
tures with thermodynamic phase transitions. In particu- 
lar, a stressed solid can be considered to be in a metastable 
state (|] and the point of global failure can be seen as a nu- 
clcation process in a first order transition near a spinodal 
Thus, the power law behavior observed experimen- 
tally in the acoustic emission before failure ^HQ] has been 
compared with the mean-field scaling expected close to a 
spinodal point ||. In analogy with spinodal nucleation 
|| , scaling behavior can only be seen when long-range in- 
teractions are present, as it is the case for elasticity, but 
should not be observable when the stress transfer function 
is short ranged. This observation is confirmed in fracture 
models with short-range elastic forces, which usually do 
not show scaling JT^,|llJ. 

Most of the theoretical investigations in this field rely 
on large scale computer simulation of lattice models where 
the elastic medium is represented by a spring (beam) net- 
work, and disorder is captured either by random dilution 
or by assigning random failure thresholds to the bonds 
jUI. The failure rule usually applied in lattice models is 
discontinuous and irreversible: when the local load exceeds 
the failure threshold of a bond, the bond is removed from 
the calculations (i.e. its elastic modulus is set to 0). Fur- 
thermore, failed bonds are never restored (no healing). 
Very recently, a novel continuous damage law has been 



introduced in lattice models |13| ] . In the framework of this 
model when the failure threshold of a bond is exceeded 
the elastic modulus of the bond is reduced by a factor a 
(0 < a < 1), furthermore, multiple failures of bonds are 
allowed. This description of damage in terms of a contin- 
uous parameter corresponds to consider the system at a 
length scale larger than the typical crack size. Computer 
simulations have revealed some remarkable features of the 
model: after some transients the system tends to a steady 
state which is macroscopically plastic, and is characterized 
by a power law distributed avalanches of breaking events. 

A very important class of models of material failure are 
the fiber bundle models (FBM) |(^,|^,^6|^,|^,|^,golj^,|2^,^J^|2^^.g; 

which have been extensively studied during the past years. 
These models consists of a set of parallel fibers having sta- 
tistically distributed strength. The sample is loaded paral- 
lel to the fibers direction, and the fibers fail if the load on 
them exceeds their threshold value. In stress controlled 
experiments, after each fiber failure the load carried by 
the broken fiber is redistributed among the intact ones. 
Among the several theoretical approaches, one simplifica- 
tion that makes the problem analytically tractable is the 
assumption of global load transfer, which means that after 
each fiber breaking the stress is equally distributed on the 
intact fibers neglecting stress enhancement in the vicinity 

of failed regions 131111111111^11 • The 

relevance of FBM is manifold: in spite of their simplicity 
these models capture the most important aspects of mate- 
rial damage and due to the analytic solutions they provide 
a deeper understanding of the fracture process. Further- 
more, they serve as a basis for more realistic damage mod- 
els having also practical importance. The very successful 
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Previous stud- 



micromechanical models of fiber reinforced composites are 
improved variants of FBM taking into account stress lo- 
calization (local load transfer) W|j20£l|,P4| , the effect of 
matrix material between fibers [ 20|J2 1| ,|2S 
possible non-linear behavior of fibers ]T 
ies of FBM addressed the macroscopic constitutive behav- 
ior, the reliability and size scaling of the global material 
strength, and the avalanches of fiber breaks preceding ul- 
timate faQure 

In this paper we generalize the FBM applying a contin- 
uous damage law for the elements in the spirit of Ref. |l3| . 
Emphasis is put on the microstructure of damage and 
its evolution with increasing load. For the case of global 
load transfer in the continuous damage model, we derive 
an exact analytic expression for the constitutive behavior 
and show that the system reaches asymptotically a steady 
state, which is macroscopically plastic. It is demonstrated 
that the continuous damage model provides a broad spec- 
trum of description of materials by varying its parameters 
and for special choices of the parameter values the model 
recovers the dry FBM and other micromechanical models 
of composites known in the literature. Next, we present a 
theoretical investigation of damaging in FBM by study- 
ing the 'dry' FBM varying the range of load transfer. In 
the case of global load transfer the model approaches the 
failure point by scaling laws, analogous to those observed 
close to a spinodal ||. However, scaling is not observed 
for local load transfer FBM 10 , ll| as it is expected for a 
spinodal instability, which can only be observed in mean- 
field theory. Increasing the range of interactions, one can 
observe that the spinodal point, defined in the global load 
transfer FBM, is approached. The evolution of damage is 
compared in the local load transfer dry and continuous 
damage FBM. Finally, we analyze the effect of quenched 
random threshold on damage localization. 

The paper is organized as follows: in section 2 we de- 
scribe the FBM and derive their constitutive behavior in 
section 3. In section 4 we discuss the local load transfer 
FMB focusing on the constitutive behaviour and on the 
cluster analysis, and in section 5 we explore the role of the 
type of disorder in the evolution of damage. Section 6 is 
devoted to discussion and conclusions. 



2 Models 

The system under consideration is composed of N fibers 
assembled in parallel on a two dimensional square lattice 
of side length L, i.e. N = L 2 . The geometrical structure 
of the model is illustrated in Fig. EL The square lattice 
corresponds to a cross section of a unidirectional fiber 
ensemble. In FBM, the fibers are considered to be lin- 
early elastic until breaking (brittle failure) with identical 
Young-modulus Ef but with random failure thresholds di, 
i = 1, . . . , N . The failure strength di of individual fibers 
is supposed to be independent identically distributed ran- 
dom variables with a cumulative probability distribution 
P(d). The fiber bundle is supposed to be loaded uniax- 
ially, and load F applied parallel to the fibers gives rise 
to a strain / of the bundle. When a fiber experiences a 



local load larger than its failure threshold the fiber fails. 
In dry FBM there is no matrix material present, which 
implies that broken fibers do not support load any more, 
and their load is redistributed to the surviving fibers. 

In the global load transfer FBM, after failure the load 
is transfered equally to all the remaining intact fibers, so 
that the load on fiber i is simply given by Fi = F/n s (F) J 
where n s (F) is the total number of surviving fibers for a 
load F. This also implies that the range of interaction be- 
tween fibers is infinite, and hence, the global load transfer 
corresponds to the mean field treatment of FBM. In the 
local load transfer FBM the load is transfered equally only 
to the surviving nearest neighbor fibers, giving rise to high 
stress concentration around failed regions (see also Fig. I]). 
We also study intermediate situations in which the load is 
transfered to a local neighborhood surrounding the failed 
fiber (i.e. a square of radius R centered on the failed fiber, 
see Fig. [I]). This model interpolates between the nearest 
neighbor local FBM and the global FBM as the range of 
interaction is increased. 

Next we generalize the model replacing the brittle fail- 
ure of fibers by a continuous damage parameter [ p"3"| . When 
the load on a fiber reaches the threshold value di the stiff- 
ness Ef of the fiber is reduced by a factor < a < 1. The 
characterization of damage by a continuous parameter cor- 
responds to describe the system on length scales larger 
than the typical crack size. This can be interpreted such 
that the smallest elements of the model are fibers and the 
continuous damage is due to cracking inside fibers. How- 
ever, the model can also be considered as the discretiza- 
tion of the system on length scales larger than the size of 
single fibers, so that one element of the model consists of 
a collection of fibers with matrix material in between. In 
this case the microscopic damage mechanism resulting in 
multiple failure of the elements is the gradual cracking of 
matrix and the breaking of fibers. 

In the following the elements of the continuous damage 
FBM will be called fibers, but we have the above two pos- 
sible interpretations in mind. Once the fiber i has failed 
its load is reduced to afi and the rest of the load (1 — a) fa 
is distributed equally among all the other fibers (global 
stress transfer) or among the neighboring fibers (local 
stress transfer). In principle, a fiber can now fail more 
than once and we define k max as the maximum number 
of failures allowed per fiber. We will first study the model 
for finite k max and eventually take the limit k max — > 00. 
It is important to note that once a fiber has failed, we can 
either keep the same failure threshold (quenched disorder) 
or chose a different one of the same distribution (annealed 
disorder), which can model microscopic rearrangements 
in the material. The failure rules of the model in the two 
cases are illustrated in Fig. [| In the following sections we 
will analyze both cases, showing that there are differences 
in the microstructure of damage between quenched and 
annealed disorder in this problem. 
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3 Constitutive laws 



constitutive equation is given by 



Here we derive the constitutive law for continuous dam- 
age FBM and show how the FBMs used in the literature 
can be recovered in particular limits. We first consider 
the case in which fibers are allowed to fail only once: the 
constitutive equation reads as 



- = f(l-P(f))+afP(f), 



(1) 



where P(f) and 1 — P(f) are the fraction of failed and 
intact fibers, respectively, and the Young-modulus Ef of 
intact fibers is taken to be unity. In Eq. ([!]) the first term 
provides the load carried by intact fibers while the second 
term is the contribution of the failed ones. Note that this 
particular case together with the parameter choice a = 
(i.e. broken fibers carry no load) corresponds to the 
dry FBM @JT|,|T^,^8| , while setting a = 0.5 in Eq. (0) 
we recover the so-called micromechanical model of fiber 
reinforced ceramic matrix composites (CMC's), which has 
been extensively studied in the literature [ p5|]26| , |27| , p8[ . In 
CMC's the physical origin of the load bearing capacity of 
failed fibers is that in the vicinity of the broken face of 
the fiber the fiber-matrix interface debonds and the stress 
builds up again in the failed fiber through the sliding fiber- 
matrix interface. 

When the fibers are allowed to fail more than once 
we have to distinguish between quenched and annealed 
disorder. 

(i) Quenched disorder: When the fibers are allowed to 
fail twice the constitutive equation can be written as 



— = f(l 
N J y 



P(f)) + af [P(f) - P(af)} + a 2 fP(af), (2) 



where [P(f) — P(af)] is the fraction of those fibers which 
failed only once, and P(af) provides the fraction of fibers 
which failed already twice. In the general case, when fibers 
are allowed to fail k max times, where k max can also go to 
infinity, the constitutive equation can be cast into the form 



F 
N 



= /(l-P(/))+ £ a*f [P(a'-V) - Pia'f)] 

+a *— /iV— - 7) (3) 



(ii) Annealed disorder: As in the previous case we con- 
sider first the case in which fibers are allowed to fail twice, 
obtaining 



§=/(l-P(/)) 



afP(f)(l-P(af)) 



(4) 



-a 2 fP(f)P(af), 



where P(f)(l — P(af)) is the fraction of fibers which failed 
only once, and P(f)P(af) is the fraction of fibers which 
failed already twice. Finally, when fibers are allowed to 
fail k max times, where k max can also go to infinity, the 



F 

N 



a i f[l-P(a i f)]l[P(aif)+ (5) 
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In Fig. |^ we show the explicit form of the constitutive 
law for quenched disorder for different values of k max in 
the case of the Weibull distribution 



F(d) = I-exp(-(d/4) m ), 



(6) 



where m is the Weibull modulus and d c denotes the char- 
acteristic strength of fibers. It is important to remark that 
the constitutive laws derived above are exact only in the 
infinite size limit (N — > oo), while fluctuations in the value 
of the failure stress F c have been observed and studied for 
finite size bundles. For this reason, we compare the theo- 
retical results with numerical simulations of bundles of size 
N = 128 2 . The agreement between simulations and theory 
turns out to be satisfactory both for quenched (Fig. |) and 
annealed disorder and reflects the fact that for global load 
sharing finite size fluctuations, for instance for F c , should 
scale as 1/ y/(N) . This behavior has to be contrasted with 
local-load sharing fiber bundles where finite size effects are 
very strong, as we will discuss in the following. 

In Fig. ||a the fibers are supposed to have a kmax resid- 
ual stiffness after having failed k max times, which gives 
rise to hardening of the material, i.e. the F/N curves 
asymptotically tend to straight lines with slope a kmam . In- 
creasing k max the hardening part of the constitutive be- 
havior is preceded by a longer and longer plastic plateau, 
and in the limiting case of k max — > oo the materials be- 
havior becomes completely plastic (see Fig. [}]). A similar 
plateau and asymptotic linear behavior has been observed 
in brittle matrix composites, where the multiple cracking 
of matrix turned to be responsible for the relatively broad 
plateau of the constitutive behavior, and the asymptotic 
linear part is due to the linear elastic behavior of fibers 
remained intact after matrix cracking p0| . 

In order to describe macroscopic cracking and global 
failure instead of hardening, the residual stiffness of the 
fibers has to be set to zero after a maximum number k* 
of allowed failures JT3]. In this case the constitutive law 
can be obtained from the general form Eqs. (|^) and (||) 
by skipping the last term corresponding to the residual 
stiffness of fibers, and by setting k max = k* in the re- 
maining part. A comparison of the constitutive laws of 
the dry and continuous damage FBM is presented in Fig. 
||b for the case of quenched disorder. Annealed disorder 
yields similar results. One can observe that the dry FBM 
constitutive law has a relatively sharp maximum, how- 
ever, the continuous damage FBM curves exhibit a plateau 
whose length increases with increasing k* . Note that the 
maximum value of F/N corresponds to the macroscopic 
strength of the material and in stress controlled experi- 
ments the plateau and the decreasing part of the curves 
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cannot be reached. However, by controlling the strain /, 
the plateau and the decreasing regime can also be realized. 
The value of the driving stress a = F/N corresponding to 
the plastic plateau is determined by the damage parame- 
ter a, while the length of the plateau is controlled by k max 
and k* . 

In Fig. U we directly compare the constitutive law for 
quenched and annealed disorder and confirm that the dif- 
ferences between the cases are very small. In particular, 
all the basic constitutive behavior are reproduced in the 
two cases. 

It is important to remark that the behavior of the dry 
FBM model (a = 0) under unloading and reloading to the 
original stress level is completely linear, since no new dam- 
age can occur during unloading-reloading sequences and 
the effect of the matrix material is completely neglected. 
This also implies that in each damage state the model is 
completely characterized by the Young modulus defined as 
the slope of the unloading curve. If the value of the dam- 
age parameter is larger than (a > 0) the behavior of the 
system under unloading and reloading is rather compli- 
cated. Due to the sliding of broken fibers with respect to 
the matrix, hysteresis loops and remaining inelastic strain 
occur (for examples see Ref. Q and references therein). 



4 Local load transfer FBM 
4.1 Constitutive behaviour 

To study the effect of stress enhancement around failed 
fibers on the damage evolution and on the macroscopic 
constitutive behavior we employ local load transfer for 
the stress redistribution after fiber failure pp| , p5H2^ , p7t . 
Since this case cannot be treated analytically, we perform 
numerical simulations in the dry FBM model: after fiber 
failure the load is redistributed on the intact nearest neigh- 
bors of the failed fiber using periodic boundary condition 
on the square lattice (see also Fig. |l|). For simplicity, in this 
case the strength of fibers di has a uniform distribution 
between and 1. The algorithm to simulate the loading 
process is as follows: (i) we impose on all the fibers the 
same load, equal to the smallest failure threshold, which 
results in breaking of the weakest element, (ii) The load 
carried by the failed fiber is redistributed on the intact 
nearest neighbors, and the load of the broken fiber is set 
to zero, (Hi) After the stress redistribution, those fibers 
whose load exceeds their failure threshold di are identi- 
fied and removed from the calculation, and the simulation 
is continued with point (ii). If the configuration obtained 
after the stress redistribution is stable, the global load is 
increased to cause the failure of one more fiber and the 
simulation is continued with point (ii). This procedure 
goes on until all fibers are broken. The applied stress just 
before global failure is considered to be the failure strength 
of the model solid. Simulations were performed with sys- 
tem sizes L = 16, 32, 64, 128. 

The constitutive behavior of the local and global load 
transfer dry FBM is compared in Fig. @. For clarity, the 



total force F (instead of stress) is presented as a func- 
tion of strain /, for several different system sizes L. Note 
that N = L 2 is chosen the same for global and local load 
sharing simulations. In the case of global load transfer the 
N — ► oo constitutive law can be obtained exactly by sub- 
stituting the cumulative probability distribution P(f) = / 
of the uniform distribution into the general form Eq. (Ty) 
and setting the damage parameter a = 0: 

F = Nf(l-P(f))=Nf(l-f), fe [0,0.5], (7) 

the strain correspondingto macroscopic failure is f c = 0.5. 
It can be seen in Fig. || that the macroscopic constitu- 
tive behavior for local load transfer always coincides with 
the global FBM solution, however, the macroscopic failure 
strength is substantially reduced in the local case giving 
rise to more brittle constitutive behavior. It is interest- 
ing to note that increasing the system size L the failure 
strength of the local FBM decreases, showing the logarith- 
mic size effect also found in the one dimensional local FBM 
[ pr| and in two dimensional fuse networks |57|, while the 
global load transfer case does not have size dependence. 

To get a deeper understanding of the behavior of the 
system as a function of the range R of load redistribution, 
we perform simulations by redistributing the load after 
fiber failure on the intact fibers in a square of side length 
2R + 1 centered on the failed fiber. The range of load re- 
distribution R is varied between 1 and (L — l)/2. Note 
that R = 1 corresponds to local load transfer on near- 
est and next-nearest neighbors, while the limiting case of 
R = (L — l)/2 recovers the 'infinite' range global load 
transfer. The comparison of the constitutive behavior in 
the local and global load transfer case is presented in 
Fig. ^. Simulations reveal that the constitutive laws ob- 
tained at different R values always fall onto the curve of 
the global load transfer case and the macroscopic failure 
strength increases with increasing range of interaction R 
approaching the strength of global FBM. For clarity, in 
Fig. ^| we indicated by vertical dashed lines the position 
of global failure at different values of R. 

To characterize the elastic response of the dry FBM 
model in a given damage state, we compute the Young 
modulus Y, defined as the slope of unloading curves as a 
function of the driving stress a = F/N: 

Y(a) = ^=l-P(f(a)) = l-f(a). (8) 

Using the constitutive law Eq. (fjj) for the global load 
transfer case, Y can be written into a closed form as a 
function of stress 

Y(a) = ±[l + VT^\. (9) 

The results on Y for global and local load transfer are 
shown in Fig. [7| where the vertical dashed lines indicate 
the position of macroscopic failure at different values of 
the redistribution range R. It can be seen that at the fail- 
ure point Y(a) has a discrete jump, the size of which de- 
creases with increasing R, but it remains finite in the limit 
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of global load transfer. Increasing R gives rise to increas- 
ing slope of Y(o) at the failure point, and in the limit 
of infinite range interaction the tangent of Y(a) becomes 
vertical at the point of failure. 



4.2 Cluster analysis 

One of the most interesting aspects of the damage mech- 
anism of disordered solids is that the breakdown is pre- 
ceded by an intensive precursor activity in the form of 
avalanches of microscopic breaking events. Under a given 
external load F a certain fraction of fibers fails immedi- 
ately. Due to the load transfer from broken to intact fibers 
this primary fiber breaking may initiate secondary break- 
ing that may also trigger a whole avalanche of breaking. If 
F is large enough the avalanche does not stop and the ma- 
terial fails catastrophically. It has been shown by analytic 
means that in the case of global load transfer the size dis- 
tribution of avalanches follows asymptotically a universal 
power law with an exponent —5/2 p 32|, however, in the 
case of local load transfer no universal behavior exists, 
and the avalanche characteristic size is bounded |ic|jll[ . 
This precursory activity can also be observed experimen- 
tally by means of the acoustic emission analysis. Acoustic 
emission measurements have revealed that for a broad va- 
riety of disordered materials the response to an increasing 
external load takes place in bursts having power law size 
distribution over a wide range J^,^,0|. 

In this section we analyze the evolution of damage in 
local load transfer FBM, comparing dry and continuous 
damage models. Instead of avalanches of fiber failures, we 
focus on the properties of clusters of broken fibers which 
are much less explored. In the following simulations we 
employ a uniform distribution for the thresholds d,. The 
load after fiber failure is redistributed on the surviving 
nearest neighbors. As the load is increased, fiber breaks 
and clusters of broken fibers are formed due to the spatial 
correlation introduced by the local load transfer. These 
clusters of broken fibers can be identified as microcracks 
formed in the plane perpendicular to the load direction. 
We monitor the damage evolution by taking snapshots of 
the clusters at different loads. In Fig. ^ the damage evolu- 
tion is shown in the dry FBM. One can observe the nucle- 
ation and gradual growth of clusters with increasing load 
F. We find that the clusters are small compared to the 
system size even before global failure (see Fig. |§|(d)), in 
accordance with the first-order transition scenario. In the 
continuous damage FBM with local load transfer the clus- 
ters of failed fibers are defined as connected sets of fibers 
having the same number of failure, taking into account 
only nearest neighbor connections. In these calculations 
we set k m ax — oo, and the simulations are stopped when 
the plastic regime is reached. 

To obtain quantitative informations on the damage 
evolution, we measure the cluster probability distribution 
n(s,F), defined as the number of clusters formed by s 
neighboring broken fibers when the applied load is F || . 
The moments (Mk(F) = J s k n(s, F)ds is the k-th mo- 
ment) of n(s, F) contain most of the information on the 



evolution of the damage. We determine n(s,F) for differ- 
ent system sizes L by averaging over the disorder. The 
total number of clusters n c = Mo as a function of the 
load is presented in Fig. |[ The increasing part of n c as 
a function of F is due to the nucleation of new microc- 
racks, and the short plateau or decreasing regime in the 
vicinity of global failure is caused by the coalescence of 
growing cracks. The inset of Fig. || demonstrates that in 
dry FBM n c obeys a simple scaling law n c — L 2 g(F/L 2 ) 
implying that the clusters are homogeneously scattered 
through the lattice. A similar scaling is observed for the 
continuous damage annealed FBM in Fig. [H)| The inset 
shows that the scaling function is independent of the dam- 
age parameter a. 

Next, we measure the average cluster size defined as 
S = MijM\ and show that it approaches a value which 
decreases with system size (Fig. [Ill a)). It can be seen 
that for a given system size the S(F) curves have two 
regimes: a slowly increasing initial regime due to the nu- 
cleation and growth of clusters, and a rapidly increasing 
part close to global failure which is caused by the coa- 
lescence of growing cracks. In dry FBM, we can simply 
rescale the data according to the law S(F, L) — s(F/L 2 ) 
and obtain a good collapse (Fig. |ll](b)). Similar results are 
obtained for the continuous damage case, but we see that 
the rescaled curves depend on a (Fig. |l2|). The larger a is, 
the smaller the clusters are, since the stress concentration 
decreases with increasing a, and the disorder gets more 
dominating. These results demonstrate that global failure 
is initiated once the crack size reaches a critical size s c af- 
ter which a crack becomes unstable. However, the average 
cluster size S does not provide a reliable estimate of s c , 
which can be obtained instead monitoring the size of the 
largest cluster S max as a function of the load. It can be 
seen in Fig. |l^ that S max reaches a value that increases 
with L, but for large L this value seems to saturate. The 
rapid increase of S max close to the failure point is due to 
the coalescence of clusters and is thus produced by a very 
small amount of fiber failures. 



5 damage localization: effect of the quenched 
disorder 

In the previous section, we analyzed the damage structure 
in the local load transfer models. In the dry FBM and in 
the case of continuous damage FBM with annealed disor- 
der we do not expect to find any non trivial damage local- 
ization for global load transfer rules, since these models 
behave effectively like in mean-field theory. On the other 
hand, quenched disorder can lead to localized structure 
for continuous damage FBM even in global load transfer 
conditions. Weak fibers are expected to fail more times 
generating an inhomogeneous damage pattern. 

In order to analyze this effect, we measure k(i), the 
number of failures at fiber i, when there is exactly one 
fiber, which has reached k max . In the bottom part of Fig. [li] 
we plot the value of k(i) and the corresponding value of 
the threshold d(i). The 'damage' k(i) shows a very irreg- 
ular pattern, which should be compared with a roughly 
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uniform structure expected for annealed disorder. In the 
upper part of the figure we display the decay of k as a 
function of d, showing how weak fibers break more often 
than strong ones. The decay is more pronounced when a 
is close to one, and becomes less important for smaller 
a. It is straightforward to obtain an analytic expression 
for k(d ) which is compared with the numerical results in 
Fig. p. 

Finally, we expect that quenched disorder should have 
an effect also on the cluster structure of local load transfer 
models. In order to confirm this point, we compare the 
number of clusters n c (F) and the average cluster size S 
for quenched and annealed disorder. The results, shown 
in Figs. [li] and [IT], indicate that damage is more localized 
when the disorder is quenched. 



6 conclusions 

We have proposed a continuous damage version of the 
FBM, that can be used to model a wide variety of con- 
stitutive behaviors. We have analyzed the development of 
damage in FBM under different conditions and compared 
the local load transfer model with different interaction 
ranges with the global load transfer model that can be 
solved exactly in the limit N — » oo. 

From the theoretical point of view, the cluster analysis 
shows the analogies between the failure in the FBM and 
nucleation in first-order phase transition. The failure point 
in global load transfer FBM plays the role of a spinodal 
point. The strain / carried by the fibers, proportional to 
the fraction of intact fibers M , close to the failure load F c 
has a diverging derivative dj jdF ~ (F c — F)^ 1 / 2 . A similar 
behavior is observed close to a spinodal instability in first- 
order phase transitions. One should note that the spinodal 
point and its associated scaling is a mean-field property, 
obtained in the limit N — ► oo, and in general is not ob- 
served for finite dimensional short-ranged models where 
nucleation occurs much before reaching the spinodal. Sim- 
ilarly, in local-load transfer (i.e. short range) FBM failure 
occurs much before the corresponding global load sharing 
instability (i.e. the spinodal) and the avalanche character- 
istic size is bounded |32|jic| , [ll| . 

We have shown that increasing the range of interaction 
the failure point is shifted towards the spinodal. A similar 
behavior is observed for instance in Ising systems when 
the range of interactions is increased |39[ . In models with 
long range stress transfer, as for instance in elastic or elec- 
tric networks, it is possible to observe the spinodal scaling 
even in finite dimensional systems • The presence of a 
spinodal instability could explain the observation of scal- 
ing properties in acoustic emission experiments ■ 

Our continuous damage FBM can reproduce a wide 
variety of elasto-plastic constitutive behaviors. A remark- 
able feature of the model is that multiple failure of the 
elements results in ductile macroscopic behavior in spite 
of the brittleness of the constituents. Similar ductile be- 
havior has been observed experimentally in fiber rein- 
forced composites made of brittle constituents [ f£2p5| , [44f . 



Experiments revealed that the mechanism of this ductil- 
ity, which is called pseudo-strain hardening, is the mul- 
tiple failure of the material [|42| , [l3p^| . Our continuous 
damage model, recovering as special cases the dry bundle 
model and micromechanical models known in the litera- 
ture, could provide a general framework for the statistical- 
micromechanical modeling of the behavior of fiber rein- 
forced composites. The fitting of experimental results in 
the framework of our model will be presented in a forth- 
coming publication. 
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Fig. 1. The geometry of the FBM. The uniaxial fiber bundle 
(left) is modeled on a square lattice (right) corresponding to 
a cross section of the specimen. The black plaquette indicates 
a broken fiber, its nearest neighborhood is shadowed, and the 
bold line shows the neighborhood for the range of interaction 
R = 2. 

a) 




f 



Fig. 2. The constitutive behavior of a single fiber of the con- 
tinuous damage model when multiple failure is allowed, (a) 
Quenched disorder: the horizontal line indicates the damage 
threshold di, which is constant in time for each fiber, (b) An- 
nealed disorder: a new threshold is extracted at random after 
each failure. 
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Fig. 3. (a) The constitutive law for the global stress transfer 
continuous damage FBM (CDFBM) for different values of k max 
for m = 2.0, a — 0.8 and quenched disorder. Plastic behavior 
is obtained in the limit of k max — > oo. (b) Comparison between 
dry and continuous damage FBM. If we allow for brittle failure 
after k* damage events, we obtain a plastic plateau followed 
by brittle failure. Symbols refer to simulations of bundles of 
size N = 128 2 and lines to the analytic calculations. 



Fig. 4. Comparison of the constitutive laws of the quenched 
and annealed case for a = 0.8, k max — 8 and k* — 8. 
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Fig. 5. The constitutive law for the global stress transfer dry 
FBM (solid line) is compared with the local stress transfer 
model. Increasing the system size the failure stress decreases. 
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Fig. 6. The constitutive law for the global stress transfer 
dry FBM (solid line) is compared with the local stress trans- 
fer FBM for different interaction ranges R. The values of 
R corresponding to the consecutive vertical dashed lines are 
1,3,5,11,15 from left to right, and the system size L = 128 
was chosen. Increasing the interaction range the failure stress 
increases approaching the value predicted by the global stress 
transfer model. 
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Fig. 8. Snapshots of the damage in the dry FBM model on 
a square lattice of size L — 128 for different values of the 
load: a) F/F c = 0.153 b) F/F c = 0.468, c) F/F c = 0.796 d) 
F/F c = 0.997. 
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Fig. 7. The Young modulus for the global stress transfer dry 
FBM is compared with the local stress transfer model for dif- 
ferent interaction ranges R. The values of R corresponding to 
the vertical dashed lines are the same as in Fig. fj. 



Fig. 9. The number of clusters as a function of the load F in 
the local stress transfer dry FBM for different system sizes. In 
the inset the rescaled plot is presented, where also the system 
L — 16 is shown (square). 




Fig. 10. The number of clusters as a function of the load F 
in the local stress transfer continuous damage annealed FBM 
for different system sizes at a = 0.5. In the inset we show the 
rescaled plot for different values of a and L. 
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Fig. 12. (a)The average cluster size as a function of the load F 
in the local stress transfer continuous damage annealed FBM 
with a = 0.5 for different system sizes L, and (b) the corre- 
sponding rescaled plot for different values of L and a. 
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Fig. 11. (a) The average cluster size as a function of the load 
F in the local stress transfer dry FBM for different system sizes 
L and the corresponding rescaled plot (b). 
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Fig. 13. The size of the largest cluster S m ax as a function of 
F/L 2 in the local stress transfer dry FBM for different system 
sizes L. 
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Fig. 14. The number of failures k of the fibers as a function 
of their threshold d (distributed uniformly) in the quenched 
global load transfer continuous FBM (top). The parameters 
used are k max = 100, N = 500 and a = 0.4, 0.8, 0.9, 0.95. The 
threshold d(i) as a function of i is compared with the 'damage' 
k(i) for a — 0.9 (bottom). 



Fig. 15. The number of failures k of the fibers as a function 
of their threshold d (distributed according to the Weibull dis- 
tribution with m = 1.5) in the quenched global load transfer 
continuous FBM: comparison between simulations and analytic 
results for a = 0.95, 0.9, 0.4. 
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Fig. 16. The number of clusters as a function of the rescaled 
load F/L in the local stress transfer continuous damage 
quenched and annealed FBM for different system sizes at 
a — 0.5. Note that the number of clusters increses faster for 
annealed disorder indicating a smaller degree localization. 




Fig. 17. The average cluster size as a function of the rescaled 
load F/L 2 in the local stress transfer continuous damage 
quenched and annealed FBM for different system sizes at 
a = 0.5. Note that the cluster size increses faster for quenched 
disorder indicating a larger degree localization. 



